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Abstract 

We  present  a  survey  of  recent  developments  In  the  theory 
of  the  Boltzmann  equation  for  a  dilute  monatomlc  gas.   These 
advances  are  characterized  by  Increased  sophistication  In 
mathematical  techniques,  Improved  understanding  of  the  signif- 
icance of  older  approximate  results,  and  a  general  trend  away 
from  ad_  hoc  and  toward  more  precise  mathematical  procedures. 
In  the  singular  limit  of  small,  mean  free  path,  the  traditional 
Hllbert  and  Chapman-Enskog  expansions  have  been  shown  to  be 
asymptotic  to  true  solutions  of  the  Boltzmann  equation,  but 
only  when  the  variables  are  appropriately  interpreted.   At 
the  opposite  extreme  of  large  mean  free  path  where  the 
behavior  Is  again  nonuniform,  the  precise  mathematical  singu- 
larities have  been  exposed.   Over  the  whole  range  of  linear 
problems  the  presence  of  a  continuous  spectrum  In  both  the 
collision  operator  and  the  streaming  operator  points  to  the 
Inadequacy  of  traditional  exponential  "normal  mode"  expan- 
sions In  both  Initial  and  boundary  value  problems.   Results 
in  all  the  reglm.es  have  been  tied  together  by  the  overall 
qualitative  understanding  given  by  a  more  comprehensive 
existence  theory  which,  for  the  first  time,  is  broad  enough 
to  encompass  the  transition  to  macroscopic  continuum  flow. 
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1 .   Introduction 

The  purpose  of  this  article  Is  to  bring  up  to  date 

1  2 
earlier  surveys  '   of  the  state  of  the  art  In  the  theory  of 

the  Boltzmann  equation.   Historically,  the  principal  methods 
of  attack  have  been  connected  with  normal .solutions  and 
transport  coefficients  (Hllbert  and  Chapman-Enskog) ,  poly- 
nomial and  other  more-or-less  arbitrary  approximations  to 
the  distribution  function.  Integral  Iterations,  and  methods 
specific  to  models  of  the  Boltzmann  equation  (e.g.  Fourier 
transform) .   Each  method  has  seen  decisive  advances  In  recent 
years.   Another  attack,  of  very  recent  appearance.  Is  direct 
numerical  computation;  this  method  Is  still  largely  explora- 
tory. 

The  significance  of  the  Hllbert  and  Chapman-Enskog 
theories  as  formal  asymptotic  expansions  of  solutions  of  the 
Boltzmann  equation  was  described  qualitatively  In  I958  and 
In  full  mathematical  detail  more  recently    (see  Sect.  3)- 
The  crucial  qualitative  point  Is  the  necessity  to  reinterpret 


1.  Harold  Grad,  "Principles  of  the  Kinetic  Theory  of  Gases"  In 
Handbuch  der  Physlk,  ed .  S.  Flugge ,  Vol. XII  (Springer 
Verlag,  1958)- 

2.  Harold  Grad,  "Theory  of  Rarefied  Gases"  In  Rarefied  Gas 
Dynamics,  ed.  P.  Devlenne  (Pergamon  Press,  19^0) . 

3-   Reference  1,  ^26. 

4.   Harold  Grad,  Phys .  Fluids  6,  l^J    (1963). 

5-   Harold  Grad,  "Asymptotic  Equivalence  of  the  Navler-Stokes 
and  Nonlinear  Boltzmann  Equations",  to  appear  In  Proc .  of 
American  Mathematical  Society  Symposium  on  Applications  of 
Partial  Differential  Equations,  New  York,  April  19^4. 
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the  variables  that  arise  naturally  In  the  Hllbert  and  Chapman- 
Enskog  expansions  as  suitably  modified  macroscopic  properties 
of  the  gas,  not  simply  the  fluid  state.  As  a  result  of  this 
generalized  Interpretation,  the  Chapman-Enskog  expansion  has 
been  found  to  be  representative  of  very  general  solutions  of 
the  Boltzmann  equation  In  the  limit  of  small  mean  free  path. 
Polynomial  and  moment  approximations,  although  to  some 

extent  forseen  by  Maxwell,  overcame  an  almost  exclusive  pre- 

6 
occupation  with  transport  coefficients  only  In  19^9   and  were 

later  extensively  developed  for  the  linear  equation  with 

7 
Maxwelllan  potential  by  Uhlenbeck  and  Wang  Chang  and  by  Mott- 

Q  Q 

Smith   (also  more  recently  by  Slrovlch  ) .   More  complicated 
and  more  accurate  (but  essentially  arbitrary)  approximations 
followed,  principally  the  half -range   '    and  the  blmodal.   ' 


6.   Harold  Grad,  Communs.  Pure  and  Appl .  Math.  2,  331  (19^9)- 

7-   G.  E.  Uhlenbeck  and  C.  S.  Wang  Chang,  "On  the  propagation 

of  sound  In  monatomlc  gases",  Univ.  of  Michigan,  Engineering 
Research  Inst.  (1952). 

8.   H.  M.  Mott -Smith,  "A  New  Approach  In  the  Kinetic  Theory  of 
Gases",  MIT,  Lincoln  Laboratory  (195^)- 

9-   L.  Slrovlch,  Phys.  Fluids  6,  2l8  (1963)- 

10.  E.  p.  Gross,  E.  A.  Jackson,  and  S.  Zlerlng,  Ann.  Phys.  ]^, 
l4l  (1957). 

11.  Lester  Lees,  "A  kinetic  theory  description  of  plane  compres- 
sible Couette  flow",  GALCIT  Hypersonic  Research  Project, 
Memo  No.  51  (1959) • 

12.  H.  M.  Mott -Smith,  Phys.  Rev.  82,  885  (I95I). 

13-   H.  Salwen,  C.  E.  Grosch,  and  S.  Zlerlng,  Phys.  Fluids  7, 
180  (1964). 
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But  only  recently  have  these  methods  been  correctly  recog- 
nized as  being  ad_  hoc  approximations  which  ignore  the 
analytical  peculiarities  of  exact  solutions.   Although  a 
judiciously  selected  low  order  approximation  may  give  a 
good  overall  description,  greater  accuracy  is  not  easily 
forthcoming  in  higher  order,  and  even  qualitative  features 
may  be  lost  if  the  precise  analytical  peculiarities 
(notably  nonuniform  behavior  and  the  continuous  spectrum) 
are  not  taken  into  account   '    (see  Sect.  5)-   In  certain 
nonlinear  flows  the  question  of  convergence  of  polynomial 

approximations  can  even  be  given  a  definite  negative  answer 

^  .    .     ,      16 
m  certain  circumstances. 

17 
Integral  iterations,  although  formulated  by  Enskog 

and  used  for  abstract  existence  theory  '  ,   have 

only  recently  received  substantial  use  in  problems  in  which 


14.  C.  Cercignani,  Ann.  Phys .  20.  219  (1962). 

15.  Harold  Weitzner,  "Steady  State  Oscillations  in  a  Gas" 
to  appear  in  Proceedings  of  the  Fourth  International 
Symposium  on  Rarefied  Gas  Dynamics,  Toronto,  1964. 

16.  L.  W.  Holway,  Jr.,  Phys.  Fluids  ]_,    911  (1964). 

17.  D.  Enskog,  Ark.  Mat.,  Astronom.  Phys.  Ser.  A  21,  1  (I928) 

18.  T.  Carleman,  Acta  math.,  Stockh.  60,  91  (1933)- 

19.  E.  Wild,  Proc.  Cambridge  Phil.  Soc.  47,  602  (I951). 

20.  D.  Morgenstern,  Proc.  Nat.  Acad.  Scl.  4o,  719  (1954). 

21.  Reference  1,  Chap.  III. 
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22 
the  method  Is  accurate,  viz.,  in  near  free  flow   ,  flow 

23 
near  a  boundary   ,  and  In  problems  where  a  good  approxlma- 

24 
tion  Is  known  as  a  first  iterate.     The  first  iterate  is 

even  surprisingly  accurate  in  problems  of  near  free  flow 

(cf.  Sect.  4). 

25  26 
As  to  models  of  the  Boltzmann  equation  -^'    a  number 

of  essentially  exact  results  have  been  obtained  by  use  of 


Fourier  and  Laplace  transform,  e.g.,  the  initial  value 

flow  £ 

1^,15 


27  28  29 

problem   ,  flow  around  a  small  object   '   ,  and  flow  in  a 


half  space 

Most  of  the  advances  listed  above  are  tied  together 
with  the  thread  of  mathematically  singular  and  nonuniform 


22.  D.  R.  Willis,  "The  Effect  of  the  Molecular  Model  on 
Solutions  to  Linearized  Couette  Flow  with  Large  Knudsen 
Number"  in  Rarefied  Gas  Dynamics,  ed .  L.  Talbot 
(Academic  Press,  19&1 ) • 

23.  S.  Zlering,  L.  Chi,  and  R.  Fante ,  "Kinetic  Theory  of  the 
Leading  Edge",  to  appear  in  Proceedings  of  Fourth  Inter- 
national Symposium  on  Rarefied  Gas  Dynamics,  Toronto, 1964 . 

24.  H.  W.  Liepmann,  R.  Narasimha,  M.  T.  Ghahine ,  Phys.  Fluids 
5,  1313  (1962). 

25.  P.  W.  Welander,  Arkiv  Fysik  7,  507  (195^)- 

26.  P.  L.  Bhatnager,  E.  P.  Gross,  and  M.  Krook,  Phys.  Rev. 
94,  511  (195^)- 

27.  L.  Sirovich  and  J.  K.  Thurber,  "Sound  Propagation  According 
to  the  Kinetic  Theory  of  Gases"  in  Rarefied  Gas  Dynamics,  ed . 

J.  A.  Laurmann  (Academic  Press,  1963) • 

28.  Harold  Grad,  "Flow  Equations  in  a  Rarefied  Atomsphere"  in 
Proc .  of  Conference  on  Aerodynamics  of  the  Upper  Atmosphere, 
Rand  Corp.  (1959) • 

29.  Marian  H.  Rose,  "Drag  on  an  Object  in  Nearly  Free 
Molecular  Flow",  to  appear  in  Phys.  Fluids,  August  1964. 
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behavior.   Although  the  physical  consequences  may  not  always 
be  great,  the  mathematical  sensitivity  which  Is  thus  revealed 
must  not  be  disregarded  In  any  attempt  to  derive  accurate 
mathematical  consequences  from  the  equations. 


2 .   General  Theory 

The  development  of  a  comprehensive  theory  of  existence 
and  qualitative  behavior  of  solutions  of  the  Boltzmann 
equation  serves  several  "practical"  purposes.   The  first  is 
to  help  one  recognize  whether  the  solution  of  a  specific 
problem  is  representative  of  more  general  cases  or  whether 
it  may  be  exceptional.   For  example,  we  recall  that,  almost 
all  the  extrapolations  that  were  drawn  from  a  very  ingenious 
explicit  solution  of  an  exact  nonlinear  Boltzmann  equation 
have  been  found  to  be  nonrepresentatlve  of  the  behavior  of 
the  typical  solution. 

A  more  basic  goal  of  a  general  theory  is  to  determine 
whether  or  not  there  exist  any  solutions  at  all.   Such  a 
failure  is  always  a  conceptual  possibility  no  matter  how 
convincing  are  the  physical  arguments  which  suggest  the 


30.   E.  Ikenberry  and  C  Truesdell,  J.  Rati.  Mech.  Anal.  5, 
1  (1956). 
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validity  of  the  equation  and  no  matter  how  plausible  are 
the  results  derived  by  nonrigorous  approximate  means.   The 
reason  is  the  inherent  sensitivity  of  mathematical  struc- 
tures to  arbitrarily  small  perturbations  which  are  undect- 
able  physically.   In  other  words,  the  validity  of  an  equation 
such  as  Boltzmann's  must  be  based  on  studies  of  internal 
mathematical  consistency  as  well  as  on  external  confrontation 
with  experiment. 

With  regard  to  the  Boltzmann  equation,  there  is  at  the 

present  time  only  the  beginning  of  a  theory  pointing  to  the 

5 
existence  of  nonlinear  solutions  in  the  large.    On  the 

other  hand,  despite  pessimistic  statements  in  the  litera- 
ti 52 
ture,   '    there  is  no  specific  evidence  to  the  contrary 

other  than  the  fact  that  it  is  too  much  to  expect  nonlinear 
equations  in  general  to  have  solutions  in  the  large  (exis- 
tence of  solutions  would  have  to  result  from  specific 
peculiarities  of  the  Boltzmann  equation).   There  are  powerful 
existence  theorems  for  spatially  homogeneous  problems  in 
which  fluid  dynamics  plays  no  role.   But  the  question  of  the 
existence  of  spatially  dependent  solutions  is  certainly  as 
difficult  as  the  corresponding  question  in  nonlinear  fluid 
dynamics,  and  this  connection  is,  in  fact,  the  most  important 


31.  D.  Morgenstern,  J.  Rati.  Mech.  Anal.  4,  533  (1955) 

32.  Ya.  Povzner,  Math.  Sbornik,  62  (I962). 
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heuristic  argument  for  a  belief  in  the  reasonable  behavior 
of  solutions  of  the  Boltzmann  equation  in  the  large. 

Turning  to  another  difficulty,  there  is  no  evidence 
at  this  time  that  even  the  linear  equation  has  any  validity 
for  potentials  of  infinite  total  cross-section.   We  would 
expect  no  physical  differences  between  an  infinite  range 
but  rapidly  decaying  intermolecular  potential  and  a  closely 
approximating  cutoff  potential;  but  the  mathematics  seems  to 
be  extraordinarily  sensitive  to  this  point  (we  must  distin- 
guish solution  of  the  Boltzmann  equation  from  solution  of 
the  integral  equation  which  arises  in  transport  theory  and 
which  is  known  to  be  insensitive  to  "small"  changes  in  the 
potential ) . 

The  first  results  on  existance  of  solutions  to  the 
Boltzmann  equation  were  obtained  by  Carleman  for  hard  sphere 

-1  o 

molecules.    He  obtained  a  strong  nonlinear  result  in  the 
large  but  only  in  a  spatially  homogeneous  problem.   Although 
the  mathematics  is  quite  difficult,  this  result  is  physically 
trivial  and,  in  particular.  Includes  no  fluid  dynamics. 
Similar  results  (spatially  homogeneous)  were  proved  much  more 

easily  for  the  mathematically  simpler  cutoff  Maxwellian 

^   ^.  ,  19,20 
potential. 

Carleman  also  gave  an  existence  proof  for  the  linear 

33 
hard  sphere  equation  Including  spatial  variation   ,  but  with 


33.   T.  Carleman,  Froblemes  Mathelnatlques  dans  la  Theorie 

Clnetique  des  Gaz   (Almqvlst  and  Wlksells,  Uppsala,  1957) 
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t/T 

a  growth  estimate  e  '       (t  is  the  mean  free  time).   This  Is  a 
catastrophically  poor  estimate  on  any  macroscopic  time  scale. 
Recently,  the  theory  of  the  linear  equation  has  been  gener- 
alized to  include  essentially  arbitrary  cutoff  repulsive 
potentials,  and  it  has  been  extended  to  show  not  only  exis- 
tence but  smoothness  and  general  good  behavior  of  solutions, 
in  particular  boundedness  for  all  time,  uniformly  in  the 
mean  free  path.    Moreover,  all  solutions  are  shown  to 
approach  equilibrium  (under  suitable  boundary  conditions).  -^ 
This  is  a  nontrivial  extension  of  Boltzmann's  H-theorem 
which,  in  its  original  formulation,  can  show  only  an  approach 
to  a  local  Maxwellian  in  a  spatially  varying  problem.   It  is 
the  fact  that  the  solution  does  not  approach  too  closely  to 
the  local  Maxwellian  which  provides  a  residual  irreversibility 
and  allows  an  eventual  fluid-like  decay  to  equilibrium. 

Based  on  this  strong  linear  theory,  a  limited  nonlinear 
theory,  for  excursions  from  linearity  which  are  not  too  large, 
has  been  obtained  for  times  of  interest  on  a  macroscopic 
fluid  dynamic  time  scale.    Thus  the  development  of  a  complete 
nonlinear  theory  would  seem  to  be  subject  to  technical  rather 
than  basic  difficulties.   It  is  noteworthy  that  solutions  of 


y\ .      Harold  Grad,  "Asymptotic  Theory  of  the  Boltzmann  Equation, 
II"  in  Rarefied  Gas  Dynamics,  ed .  J.  A.  Laurmann,  Vol.  I 
(Academic  Press,  19fc'3)  • 

35-  Harold  Grad,  "On  Boltzmann's  H-theorem",  to  be  published 
in  Proceedings  of  SIAM  von  Karman  Symposium,  Washington, 
May  1964. 
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the  linear  Boltzmann  equation  approximate  the  true  nonlinear 
solutions  better  than  might  be  expected  on  the  basis  of 
plausible  estimates. 

For  certain  altered  versions  of  the  Boltzmann  equation 

chosen  for  mathematical  convenience,  strong  nonlinear 

31  32 
existence  theorems  In  the  large  have  been  proves. 

But  these  equations  contain  arbitrary  averaging  or  smoothing 
operators,  and  the  relation  to  an  actual  gas  is  not  estab- 
lished.  These  replacements  for  the  Boltzmann  equation 
should  not  be  confused  with  so-called  "model"  equations 
which  do  retain  a  heuristic  physical  meaning  but  do  not 
appreciably  simplify  nonlinear  existence  theory. 


3-   Small  Mean  Free  Path 

On  the  basis  of  physical  evidence  we  expect  the 
detailed  description  of  a  gas  by  means  of  a  molecular 
distribution  function  to  be  adequately  approximated  by  an 
apparently  cruder  macroscopic  fluid  description  when  the 
mean  free  path  is  small.   Mathematically  this  must  corres- 
pond to  a  very  singular  limit.   The  appropriate  mathematical 

36 
formulation,  introduced  by  Hllbert   ,  is  to  place  a  large 


36.  D.  Hllbert,  Grundz-Qge  einer  allgemelnen  Theorie  der 

linearen  Integralgleichungen  (Wien,  J.  B.  Teubner,  1924 ) 
(also  Chelsea  Publishing  Co. ,  New  York,  1953)- 
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factor  l/e  (e  Is  a  scale  factor  for  the  mean  free  path  or 
mean  collision  time)  In  front  of  the  quadratic  collision 
term  in  the  Boltzmann  equation 


|I+  |.  Vf  =  I  Q(f,f)  .  (3.1) 


Despite  the  fact  that  e  =  0  is  singular,  Hilbert  looks  for 
solutions  in  the  form 

f  =  f°  +  ef^  +  e^f^  +  . . .  .  (5.2) 

Inserting  the  power  series  and  equating  like  terms,  we  find 

n 
that  the  successive  terms  f  are  determined  uniquely  by  an 

algorithm  (successiA/ely  solving  integral  and  partial  differ- 
ential equations)  which  requires  specification  of  fluid 
initial  data  only.   We  can  interpret  this  result  as  a  state- 
ment that  the  expansion  (3-2)  is  valid  for  a  special  class 
of  solutions  which  has  the  property  that  the  distribution 
function  f  at  a  given  instant  is  uniquely  determined  by  the 
instantaneous  fluid  state  (as  a  functional,  not  Just  a  point 
function  since  the  f  depend  on  derivatives  of  the  fluid 
variables).   This  observation  leads  naturally  to  the  Chapman- 
Enskog  theory.   If  we  evaluate  the  special  Hilbert  f  in  its 
dependence  on  the  fluid  state,  we  can  compute  from  it  the 
stress  tensor  and  heat  flow  vector  as  moments  of  f ;  they  will 
also  be  expressed  as  power  series  in  e  involving  successively 
higher  derivatives  of  the  fluid  variables.   The  Boltzmann 

-  13  - 


equation  can  now  be  bypassed,  and  the  time  evolution  of 
the  fluid  state  can  be  obtained  from  the  equations  of 
conservation  of  mass,  momentum,  and  energy  which  are  self- 
contained  In  the  fluid  variables  (to  any  given  order  of  e). 

At  this  state  both  the  Hllbert  and  Chapman-Enskog 
expansions  are  only  formal;  the  exact  relation  to  a  hypo- 
thetical special  class  of  solutions  or  possibly  to  general 
solutions  of  the  Boltzmann  equation  remains  to  be  discovered. 
Also,  whereas  the  Hllbert  algorithm  can  be  carried  out  step 
by  step  by  solving  equations  with  an  established  theory,  the 
more  desirable  Chapman-Enskog  expansion  requires  the  solution 
of  nonstandard  differential  equations  (at  least  beyond  the 
Navler-Stokes  level). 

To  fix  the  Ideas  let  us  now  return  to  direct  solution 
of  the  Boltzmann  equation  without  recourse  to  power  series 
and  look  at  a  problem  with  fixed  molecular  boundary  condi- 
tions (specular,  diffuse,  or  any  other)  and  a  fixed  Initial 
function  f(0)  which  Is  Independent  of  e .   We  suppose  that 
the  equation  (3-1)  has  a  solution  f(^,x,t;  e)  and  pose  the 
question,  what  Is  the  relation  between  the  exact  solution 
and  the  Hllbert  or  Chapman-Enskog  expansions?   We  see 
Immediately  that  the  expansion  (3-2)  cannot  solve  this  problem. 
For  one,  a  Hllbert  function  (3-2)  Is  explicitly  and  uniquely 
determined  In  Its  dependence  on  molecular  velocity  ^  whereas 
f(0)  should  be  arbitrary.   Even  If  we  attempt  to  restrict  the 
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validity  of  (3.2)  to  special  Initial  values  f(0)  of  the 
Hllbert  type,  we  are  blocked  by  the  fact  that  a  Hilbert 
function  (3.2)  is  necessarily  a  function  of  e  whereas  f(0) 
has  been  chosen  independent  of  e.   The  resolution  of  this 
problem  has  been  obtained  in  recent  years  by  first  proving 
the  existence  of  the  exact  solution  f(^,x,t;  e)  (for  the 
linear   and  slightly  nonlinear-'^  equations,  see  Sect.  2); 
and  second  by  proving  that  with  appropriate  initial  date 
the  expansion  (3.2)  is  asymptotic  to  the  exact  solution. 
Specifically,  the  local  Maxwellian  f  with  fluid  parameters 
taken  from  the  solution  of  the  inviscid  Euler  equations  of 
fluid  dynamics  yields  the  result 


11m  f(l,x,t;  e)  =  f°(|,x,t)  (3.3) 

e— >  o 

at  any  fixed  positive  time  t  >  0.   Similarly,  the  Chapman- 

Enskog  second  approximation  f  +  ef  with  fluid  variation 

taken  from  the  Navier-Stokes  equations  satisfies 


lim  -  I  f  -  (f°  +  ef^)  I  =  0  (3-^) 

'       £->  O  ^ 

uniformly  in  t  >  t  >  0  (t  bounded  away  from  zero)  provided 
that  the  Navier-Stokes  equations  are  given  proper  initial 
values  --  not  the  actual  initial  fluid  state  as  specified  by 
f(0).   The  asymptotic  series  (3.2)  takes  over  only  after  an 
initial  transient  (duration  comparable  to  e).   The  lasting 
effect  of  this  transient  on  the  fluid  state  is  taken  into 
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account  by  providing  the  Chapman -Enskog  theory  with  an 
asymptotic  initial  fluid  state  which  is  derived  from  the 
entire  initial  state  f (O)  in  a  complex  way  and  is  dependent 
on  e  even  though  the  actual  Initial  fluid  state  is  not. 

This  theory  is  mathematically  complete  for  the  linear 
and  slightly  nonlinear  initial  value  problem;  an  algorithm 

has  been  given  for  the  correct  asymptotic  fluid  initial 

4 
values  for  the  general  nonlinear  problem.    By  a  plausible 

heuristic  argument   one  can  show  in  the  general  nonlinear 
problem  that  the  expansion  (3-2)  should  be  asymptotic  to  a 
true  solution  provided  that  one  avoids  not  only  the  vicinity 
of  t  =  Oj,  but  also  the  boundary  of  the  domain,  and  any 
shocks  that  may  appear  (initial,  boundary,  and  shock 
layers).   Because  these  expansions  are  only  asymptotic,  we 
are  faced  with  the  question  of  matching  the  expansions 
across  the  layers.   Only  the  problem  of  the  initial  layer 
described  above  has  been  solved  (initial  slip).   The  boundary 
layer  problem  (boundary  slip)  is  to  find  proper  asymptotic 
boundary  values  for  the  Chapman-Enskog  equations.   Any 
special  solution  of  the  Boltzmann  equation  (e.g.  heat  flow 
between  parallel  plates)  would  yield  as  a  byproduct  the 
correct  asymptotic  boundary  value  for  the  special  geometry. 
But  only  ad  hoc  approximations  to  solutions  of  boundary 
value  problems  have  thus  far  been  obtained  (with  the  excep- 
tion of  a  few  solutions  using  a  m.odel  of  the  Boltzmann 
equation  '    -^) ,    see  Sect.  5- 
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The  problem   of  shock  slip   is  to  find  the  correct 
jump  conditions  to  replace  the  Hugoniot  relations  which 
are  valid  only  to  lowest  order  as  appropriate  to  the  invis- 
cld  Euler  equations.   It  is  a  misleading  fact  that,  whereas 
the  Euler  equations  admit  only  a  discontinuity  across  a 
shock,  the  Navier-Stokes  equations  yield  a  continuous  transi- 
tion.  The  Navier-Stokes  solution,  even  though  it  is  smooth, 
does  not  correctly  approximate  the  Boltzmann  solution 
uniformly  across  a  shock  wave.   A  treatment  of  the  shock 

slip  problem  using  the  Navier-Stokes  equations  in  the  shock 

"57 
interior  '  is  properly  interpreted  as  only  a  rough  approxi- 
mation in  general,  but  one  which  should  be  expected  to  give 
the  correct  result  for  weak  shocks.   The  strictly  correct 
usage  of  the  Navier-Stokes  equations  as  they  relate  to  the 
limit  of  a  Boltzmann  solution  for  small  mean  free  path  is 
to  solve  using  the  Navier-Stokes  equations  on  both  sides  of 
a  discontinuous  shock  across  which  modified  Hugoniot  conditions 
(to  order  e)  are  applied. 

A  numerical  analysis  of  the  classical  Riemann  problem 

38 
of  a  ruptured  diaphragm   shows  qualitative  agreement  with 

this  theory.   After  an  initial  definitely  n on -Maxwell Ian 


37.  Ronald  P.  Probstein  and  Y.  S.  Pan,  "Shock  Structure  and  the 
Leading  Edge  Problem",  in  Rarefied  Gas  Dynamics,  Vol.  II, 
ed .  J.  A.  Laurmann  (Academic  Press,  19^3) • 

38.  C.  K.  Chu,  "A  Kinetic -Theoretic  Description  of  the  Forma- 
tion of  a  Shock  Wave",  to  appear  in  Phys.  Fluids. 
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transient  J  there  emerge  a  moderately  non -Maxwell Ian  shock 
layer  and  a  closely  Maxwelllan  entropy  layer.   Quantitative 
agreement  is  not  yet  feasible. 


^ .   Large  Mean  Free  Path 

This  limiting  case  (large  e)  also  exhibits  nonuniform- 
ities  but  not  so  radically  as  in  the  continuum  limit  (small 
e)  since  the  description  of  the  state  of  the  gas  in  the 
limit  of  free  flow  is  by  a  distribution  function  just  as 
for  finite  mean  free  path.   The  nonuniform  behavior  arises 
from  two  sources:   first  the  universal  presence  of  a  residual 
class  of  particles  which  is  dominated  by  collisions  no  matter 
how  large  is  the  mean  free  path;  second,  in  the  case  of  a 
molecule  of  infinite  total  cross-section  only,  the  noncon- 

vergence  of  the  collision  integral  for  a  discontinuous 

39 

distribution  function  of  the  type  occurring  in  free  flow.  ^ 

To  be  specific,  consider  a  flow  bounded  by  two  parallel 
plates  (e.g.  heat  flow  or  shear  flow).   A  particle  oriented 
in  a  grazing  direction  to  a  plate  will  travel  m.any  mean  free 
paths  between  plates  even  if  the  normal  spacing  is  much  less 


39.   We  may  note  that  although  the  infinite  cross-section 
causes  no  difficulty  in  the  formal  evaluation  of  the 
Hllbert  and  Chapman -En  skog  algorithms  for  sm.all  e,  the 
proofs  seem  to  be  very  difficult. 
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than  a  mean  free  path.   Thus  the  distribution  function  will 
differ  from  Its  free  flow  value  by  a  perturbation  of  order 
l/e  (e  is  the  Knudsen  number  and  is  large)  except  for  a 
small  region  in  velocity  space.   This  nonuniform  behavior 

is  exhibited  In  an  expansion  of  macroscopic  variables  (e.g. 

22 
heat  flow)  which  takes  the  form 

q^  +  (|  log  e)q^  +  ...  (4.1) 

where  q  is  the  free  flow  value  of  the  heat  flow.   In  this 
^o 

problem  the  lowest  order  perturbation  of  the  molecular  dis- 
tribution from  free  flow  is  correctly  obtained  uniformly  in 
velocity  by  a  first  Integral  iteration.   In  other  words,  the 
correct  distribution  Is  obtained  by  a  single  iteration 

applied  to  the  free  flow  solution.   It  is  easy  to  verify 

-2 
that  the  singularity  enters  as  e   log  e  in  a  cylindrical 

-3 
problem  and  e   log  e   in  a  domain  of  finite  extent  because 

the  group  of  exceptional  particles  is  of  successively  smaller 
dimension. 

An  entirely  different  type  of  near  free  flow  problem  is 
typified  by  flow  around  a  small  object  (or  flow  through  a 
small  orifice).   For  deflniteness  let  us  take  a  body  of  finite 
size  and  allow  the  mean  free  path  to  grow.   At  any  fixed  dis- 
tance from  the  body  the  flow  will  approach  the  classical  free 
flow,  separately  given  In  the  cone  subtended  by  the  body  and 
its  exterior.   But  this  is  evidently  not  the  case  one  mean 
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free  path  or  farther  from  the  body;  the  free  flow  solution 

p  Q  Ur\ 

is  not  uniformly  valid  out  to  Infinity.   '    To  find  the 
perturbation  uniformly  for  all  distances  Is  a  difficult 
problem.   Further,  the  leading  term  in  l/e  for  the  distribu- 
tion function  is  not  correctly  given  by  a  single  iteration  as 

41 
in  the  case  of  the  parallel  plates.   However,  it  appears 

that  the  perturbed  distribution  at  the  body  is  correctly 
given  by  one  iteration  (successive  iterations  would  presum- 
ably improve  the  accuracy  away  from  the  body).   Thus  the 
drag,  which  requires  knowledge  of  f   only  at  the  body,  can 
presumably  be  correctly  found  by  a  single  Iteration.   The 

entire  perturbed  distribution  has  been  evaluated  for  a  model 

29 
of  the  Boltzmann  equation  and  the  drag  computed.     The  first 

Integral  iteration  has  also  been  computed,  but  for  a  different 

42 
model   ,  so  the  results  cannot  yet  be  compared. 

There  are  two  sources  of  nonunif ormity  in  the  expansion 

with  respect  to  e  of  the  flow  around  an  object.   In  addition 

to  the  effect  of  the  slow  particles  as  in  the  problem  of 

parallel  plates  (or  more  precisely  concentric  spheres),  there 


40.  Y.  Pao,  "A  Uniformly  Valid  Asymptotic  Theory  of  Linearized 
Rarefied  Gas  Flows  Under  Nearly  Free  Molecular  Conditions" 
to  appear  in  Proceedings  of  Fourth  International  Symposium 
on  Rarefied  Gas  Dynamics,  Toronto,  1964. 

41.  Conversation  with  D.  R.  Willis  and  Y.  Pao. 

42.  G.  J.  Maslach,  D.  R.  Willis,  S.  Tang,  and  D.  Ko. ,  "Recent 
Extension  of  Nearly  Free  Molecular  Flow  Experiment  and 
Theory",  to  appear  in  Proceedings  of  Fourth  International 
Symposium  on  Rarefied  Gas  Dynamics,  Toronto,  1964. 
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is  also  the  effect  of  distant  collisions  on  particles  of  all 
apeeds.   Both  are  higher  order  effects  which  have  not  yet 
been  computed. 

The  nonunlf ormlty  Introduced  by  a  power  law  potential 
of  Infinite  total  cross-section  has  been  evaluated  formally, 
but  so  far  without  a  proof  of  convergence;  the  singularity 
in  a  parallel  plate  geometry  turns  out  to  be  a  fractional 
power  of  e. 


5-   Influence  of  the  Spectrum 

Essentially  all  the  difficulty  of  the  linear  Boltzmann 
equation 


II  +  ^^f  +  L(f)  =  0  (5.1) 


arises  from  the  interplay  between  the  streaming  term  i  ■  Vf 
and  the  collision  term  L(f).   They  can  be  considered  individ- 
ually by  studying  free  flow 


11+  ^-Vr  =   0  (5.2) 


4^.   J.  J.  Smolderin,  "Near  Free  Molecular  Behavior  of  Gases 
with  Infinite  Collision  Cross-section",  to  appear  in 
Proceedings  of  Fourth  International  Symposium  on 
Rarefied  Gas  Dynamics,  Toronto,  1964. 
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and  the  spatially  homogeneous  problem 


^ 


+  L(f)  =  0  .  (5-3) 


The  theory  of  Llouvllle's  equation  (5-2)  Is  dominated  by 
the  fact  that  the  spectrum  of  the  streaming  operator  |-V 
Is  continuous  and  comprises  the  entire  Imaginary  axis.   As 

a  result,  there  Is  an  approach  to  equilibrium  for  any  Initial 

44 
distribution  f(0)   ,  but  It  Is  not  universal,  and  the  rate  of 

45 
decay  to  equilibrium  may  be  exponential   or  It  may  not  and  it 

44 
will  depend  on  the  specific  initial  function  f(0). 

The  collision  operator  L  has  also  been  found  to  have  a 

continuum  as  part  of  its  spectrum  for  almost  all  cutoff 

34 
Intermolecular  potentials.     To  be  precise,  L  can  be  decom- 
posed 

L(f)  =  vf  -  K(f)  (5.^) 

where   v(^)   is  the  collision  frequency  of  a  molecule  of 
speed  i    and  K  is  nonsingular  and  compact.   The  continuous 
spectrum  of  L  consists  of  the  set  of  values  taken  by  v(^)  as 
I  ranges  over  all  Its  values.   For  a  hard  potential  v  -^  ^ 
for  large  ^ -,    therefore,  the  continuum  extends  from  some  finite 
value  to  infinity.   For  a  soft  potential  v  — >  0  for  large  i, 
and  the  continuum  extends  from  a  finite  value  to  the  origin. 


44.  Reference  1,  ^12;  and  Harold  Grad,  Communs.  Pure  and 
Appl.  Math.  XIV,  323  (l96l),  ^   5- 

45.  H.  L.  Frisch,  Phys.  Rev.  IO9,  22  (1958). 
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There  may  also  be  point  eigenvalues  and  there  is  always  a 
multiple  eigenvalue  at  the  origin  corresponding  to  conserva- 
tion of  mass,  momentum,  and  energy.   Qualitatively,  the 
relaxation  to  equilibrium  in  (5-5)  will  not  be  a  simple 
superposition  of  normal  modes. 

The  nature  of  the  spectrum  when  streaming  and  collisions 
are  coupled  is  more  complex.   For  simplicity  consider  an 
initial  value  problem  for  a  single  wave  number  k, 

11+  i(l-k)f  +  L(f)  =  0  .  (5.5) 

Now  the  relevant  operator  is  i(|-k)  +  L  and  the  continuum  is 

46 
the  set  of  complex  values  taken  by  i(4-k)  +  v(^).     Also, 

the  point  eigenvalues  of  L  at  the  origin  are  perturbed  to 

46 
neighboring  points  in  the  right  half  plane.     If  the  inter- 
molecular  potential  is  hard  and  k  is  small,  the  asymptotic 
behavior  after  a  long  time  will  be  dominated  by  these  point 
eigenvalues  and  will  approximate  the  Navier-Stokes  asymptotic 
behavior.   If  the  potential  is  soft  or  if  k  is  large,  the 
continuum  will  dominate  the  ultimate  behavior  which  will  bear 
no  resemblance  to  fluid  dynamics.   However,  if  k  is  small, 
even  with  the  soft  potential  (continuum  extending  to  the  origin) 
there  will  be  an  intermediate  time  during  which  the  point 
spectrum  dominates  (fluid-like  behavior),  although  the  ultimate 


46.   This  is  a  plausible  mathematical  conjecture  but  is  not  a 
consequence  of  known  theorems. 
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complicated  behavior  will  be  dominated  by  the  continuum. 

The  situation  Is  not  quite  the  same  In  a  boundary  value 
problem.   For  example,  In  a  half  space  we  distinguish  time- 
Independence 

(^  V)f  +  L(f)  -  0  (5-6) 

(e.g.  heat  flow  or  shear  flow)  from  the  case  of  a  vibrating 
wall 

Icof  +  (^- V)f  +  L(f)  =  0  (5-7) 

(sound  propagation  at  given  cd)  .   In  these  cases  the  continuum 

will  consist  of  the  values  taken  by  v{^)/i      or   (v+loj)/^ 

46 
respectively.     In  either  case  the  continuum  extends  to  the 

origin  (v/^  — >  0  for  large  i)   with  the  single  exception  of 
the  rigid  sphere  molecule.   In  the  steady  problem  (5-6),  the 
degenerate  point  eigenvalue  at  the  origin  dominates  far  away 
from  the  wall  and  gives  a  linear  profile  of  temperature  and 
velocity.   But  the  continuum  which  extends  to  the  origin 
implies  that  the  rate  at  which  the  solution  approaches  this 
linear  profile  is  not  a  simple  exponential.   With  the  oscil- 
lating wall  (5.7),  the  continuum  always  dominates  sufficiently 
far  from  the  wall,  and  the  asymptotic  behavior  has  no  relation 
to  the  fluid  (Navier -Stokes)  solution  of  the  same  problem. 
More  important,  as  a  consequence  of  the  continuum,  the  ultimate 
decay  will  not  be  a  property  of  the  gas  alone  (as  it  is  for  a 
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point  eigenvalue)  but  will  depend  on  details  of  the  boundary- 
condition  at  the  wall.   On  the  other  hand.  If  k  Is  small, 
there  will  be  an  Intermediate  region  where  the  displaced  point 
spectrum  (only  slightly  damped)  dominates  and  Is  approximately 
fluid  dynamic . 

These  boundary  value  problems  have  been  successfully 
solved  for  a  model  of  the  Boltzmann  equation  (for  shear  flow 
by  Cerclgnanl   and  for  sound  dispersion  by  Weitzner  ) and  the 
qualitative  features  just  described  (which  should  remain  true 
for  the  actual  Boltzmann  equation)  have  been  found. 

Perhaps  the  m.ost  significant  feature  of  the  appearance  of 
a  continuous  spectrum  is  that  there  is  no  universal  asymptotic 
decay  (either  in  time  or  in  space),  but  the  rate  of  decay 

depends  on  the  smoothness  of  the  data  (initial  or  boundary 

hi 
respectively   ) .   Quantitative  solution  of  the  actual 

Boltzmann  equation  in  boundary  value  problems  is  a  difficult 
matter.   But  even  this  qualitative  conclusion  may  turn  out  to 
be  verifiable  experimentally.   In  an  experiment  with  an  oscil- 
lating crystal,  if  co  is  taken  to  be  comparable  to  or  higher 
than  the  mean  collision  frequency,  one  should  observe  an  asymp- 
totic decrement  which  is  not  reproducible  if  the  preparation  or 
composition  of  the  crystal  surface  is  altered.   This  would 
qualitatively  contradict  not  only  the  Navler-Stokes  solution 


47.   In  a  neutral  gas  see  Ref.  44)  for  the  equivalent  problem 
in  a  plasma  ("Landau"  damping),  see  Harold  Weitzner,  Phys 
Fluids  6,  1123  (1963). 
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but  any  solution  arising  from  a  single  ultimately  dominant 
normal  mode . 

With  regard  to  potentials  of  Infinite  total  cross- 
section,  only  one  result  Is  known,  viz.,  that  there  Is  a 

pure  point  spectrum  for  Maxwelllan  molecumes  In  the  spatially 

7 
homogeneous  problem.    What  the  spectrum  Is  In  a  spatially 

dependent  problem  for  the  Maxwelllan  potential  Is  not  certain, 

and  the  spectrum  for  any  other  Infinite  range  potential  even 

In  the  spatially  homogeneous  case  is  Impossible  to  guess 

because  the  cutoff  Maxwelllan  is  borderline  and  separates 

the  case  of  a  continuum  which  extends  to  infinity  from  the 

case  of  a  continuum  which  extends  to  the  origin. 


6.   Conclusions 

Recent  developments  in  the  theory  of  the  Boltzmann 
equation  indicate  a  trend  away  from  ad  hoc  and  heuristic 
approximation  methods  and  towards  more  precise  treatments 
which  pay  more  attention  to  the  mathematical  idiosyncracies 
of  individual  problems.   There  are  a  number  of  directions  in 
which  one  can  hope  to  see  further  advances  of  this  type 
within  the  near  future. 

First  we  mention  the  influence  of  a  continuous  spectrum 
which,  although  possibly  even  qualitatively  observable 
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experimentally  (cf.  Sect. 5)  requires  quantitative  exploita- 
tion.  Up  to  now,  this  has  only  been  done  with  models  of  the 
Boltzmann  equation  which  do  not  permit  quantitative  experi- 
mental comparisons. 

A  rigorous  theory  of  the  Boltzmann  equation,  which  has 
thus  far  been  successful  only  for  cutoff  potentials,  should 
be  extended  to  cases  of  Infinite  total  cross-section.   The 
principal  purpose  Is  mathematical,  viz.,  to  find  out  whether 
such  a  very  singular  form  of  the  Boltzmann  equation  Is  a 
legitimate  mathematical  enterprise.   For  the  same  reason  there 
Is  need  for  the  development  of  strong  nonlinear  existence 
theorems  even  for  cutoff  potentials. 

The  asymptotic  interpretation  of  the  Hilbert  and  Chapman- 
Enskog  theories  requires  completion  to  Incorporate  proper 
shock  slip  and  boundary  slip  conditions.   Although  these 
problems  are  more  difficult  than  the  already  solved  initial 
slip  problem,  they  should  not  be  totally  inaccessible. 

At  the  opposite  extreme,  almost  free  flow,  the  accuracy 
of  the  first  Integral  iteration  suggests  that  extension  of 
precise  results  from  model  equations  to  the  actual  Boltzmann 
equation  will  be  quickly  forthcoming. 
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